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The paper continues the analysis of vacuum Rabi oscillations we started in Part I [Phys. Rev. 
A 79, 033836 (2009)]. Here we concentrate on experimental consequences for cavity QED of two 
different classes of representations of harmonic oscillator Lie algebras. The zero-temperature mas- 
ter equation, derived in Part I for irreducible representations of the algebra, is reformulated in a 
reducible representation that models electromagnetic fields by a gas of harmonic oscillator wave 
packets. The representation is known to introduce automatic regularizations that in irreducible 
representations would have to be justified by ad hoc arguments. Predictions based on this represen- 
tation are characterized in thermodynamic limit by a single parameter responsible for collapses 
and revivals of Rabi oscillations in exact vacuum. Collapses and revivals disappear in the limit 
<; — > CO. Observation of a finite c would mean that cavity quantum fields are described by a non- 
Wightmanian theory, where vacuum states are zero-temperature Bose-Einstein condensates of a 
A'^-particle bosonic oscillator gas and, thus, are non-unique. The data collected in the experiment 
of Brune et al. [Phys. Rev. Lett. 76, 1800 (1996)] are consistent with any ? > 400. 

PACS numbers: 42.50.Lc, 42.50.Dv, 32.80.Ee, 32.80.Qk 



I. INTRODUCTION 

The first part of this work [l| was devoted to the problem of understanding tiie structure and origin of decoherences 
that occur if a single atom propagates through an initially empty optical cavity. We compared theoretical predictions 
with the most precise available data and showed that some of the controversial issues discussed earlier in the 
literature can be resolved if one correctly takes into account the Gaussian structure of the cavity mode and the fact 
that the cavity was open. 

As opposed to standard cavity QED approaches, modeling energy losses by jumps between bare energy eigenstates 
,0], we based our analysis on Markovian master equations whose Davies operators describe jumps between dressed 
states. Our approach is consistent with the general formahsm of Davies [1,11], or the more recent works of Scala et 
al. H, Q , where losses due to interaction with environment lead to transitions between eigenstates of the system 
Hamiltonian (guaranteeing stationarity of the asymptotic state at T > K). What is important, the choice of dressed 
states simplifies computations and naturally incorporates long-wavelength transitions within a single dressed-state 
manifold, an effect expected in open cavities. In order to derive the model from a microscopic level, we assumed the 
system-reservoir interaction of the form (a(a + a^) + fia^ a) ® B, where a, (3 are parameters. 

The part proportional to j3 (let us refer to it as the Alicki interaction term [13]) is responsible for transitions 
within the same dressed-state manifold. It simultaneously makes the coupling between the system and the reservoir 
more sensitive to the photon number. Perhaps, it is the latter property of the interaction that leads to apparent 
underestimates of the cavity quality factor if Q is measured with relatively strong fields, while the actual measurements 
of Rabi oscillations are performed in almost exact vacuum. The problem of Q is one of the issues that require further 
experimental and theoretical studies. 

Still, the list of open questions is longer, and some of them touch the very fundamentals of quantum field theory. 
For example, it is known that different physical systems in general correspond to different representations of Lie 
algebras. Fields are quantized by means of harmonic-oscillator Lie algebras, but can the cavity QED data tell us 
something about their representations? It turns out that there exists a class of physically motivated representations 
whose predictions are characterized, in certain thermodynamic limit, by a single parameter which influences vacuum 
Rabi oscillations. 

Determination of <; is, in principle, within the reach of cavity QED experiments. The first estimates on <7, > 200, 
were given in [ll| , but the approach to decoherence was not in that paper based on systematically derived master 
equations. As such, it was not reliable, a fact that motivated the research project whose partial results were reported 
in [l|, and now completed in the present paper. Basing the analysis of decoherence on the results from we will 



2 



show that the data from the experiment of Brune et al. [2*1 are consistent with any <^ > 400. We will also show that 
a finite value of <r implies collapses and revivals of Rabi oscillation even in exact vacuum. The first revival time is 
tr = {(^ + A/c(<r — l))TRabi, whcrc Tptabi IS the period of Rabi oscillation. 

Observation of the revival in exact vacuum could be a proof of a non-Wightmanian nature of cavity QED. Conceptual 
consequences of such a finding might be enormous but — paradoxically — implications for agreement between standard 
theory and experiment could be smaller from what one might expect at a first glance. The reason is the correspondence 
principle stating that the weak law of large numbers, N —> oo with Z — const, maps theories based on our reducible 
representations into regularized forms of those based on irreducible representations. The limiting forms are already 
regularized, that is, the automatic cut-offs occur in exactly those places where in standard approaches one puts them 
by hand. This is a strong argument in favor of field quantization in terms of reducible representations of harmonic 
oscihator Lie algebras. 

The paper is organized as follows. In Sec. II we discuss the physical background of the reducible representation in 
question. In Sec. Ill we show how to decompose the reducible representation into blocks that allow us to perform 
calculations, in each block separately, by means of the methods known from the standard formalism. In Sec. IV we 
describe in detail the block structure of dressed states. We show that the number of dressed states in a given block 
is s + 2, where s is a parameter that characterizes the block. In Sec. V we generalize to the reducible representation 
the derivation of an appropriate master equation; the strategy is the same as in [1], only the representation of the 
Lie algebra is different. As opposed to fH, we concentrate on zero-temperature master equation. Therefore, we had 
to supplement this section be a technical Appendix on T = irreducible-representation solutions (given at the end 
of the paper). An extension to T > would be immediate, but it would simultaneously introduce a number of 
irrelevant technical details, thus obscuring the main message of the paper. In Sec. VI we discuss energy losses and 
compare the reducible treatment with the one based on general irreducible representations. The goal of this section 
is to give an alternative proof of necessity of renormalizing the decay parameters occurring in the master equations, 
before one compares theory with experiment. Sec. VII contains the main result of the paper: experiments can, in 
principle, discriminate between irreducible and reducible representations, and we show how to estimate the relevant 
experimental parameters. 



II. PRELIMINARIES 

In order to understand the wider context of the issue, let us return to the Hamiltonian we employed in [l|, 

H — (73 + o^a^a + (7((T_a^ + (T-l_a)^ 

+ {a{a + a^)+(3a^a)B + HR. (1) 

CTa, a± = (cti ± z(T2)/2 are the Pauli matrices, a is the usual harmonic-oscillator annihilation operator, the coupling 
parameter q is, for simplicity, assumed to be real, and B and Hji are operators whose explicit form is irrelevant since 
they correspond to the reservoir. It is known, at least since the work of Tavis and Cummings jl^], that ^ is the 
simplest case of a more abstract Hamiltonian 



H 



h(^ujJ3 + Loa^a + q[J-a) + J+a)^ 

+ {a{a + a'^)+ l3a'^a)B + HR, (2) 

where J3, J± are elements of the Lie algebra su(2), 

[J3,^±] = ±J±, (3) 

[J_,J+] = 2J3. (4) 

Replacing a single two-level atom by a system of several two-level atoms, one finds that J3, J± are given by a 
higher-spin, reducible representation of su(2). 

From the point of view of a Lie-algebraic purist, the Hamiltonian ([2]) mixes abstract elements J3, J± of su(2) with 

a concrete representation of another Lie algebra. Indeed, the operators a_ = a, a+ = a\ 03 — a^a, and — \ are a 
representation of the one-dimensional harmonic-oscillator Lie algebra [l3| . ho(l), 

[a_,a+] = ao, (5) 

[a±,ao] = [03, ao] = 0, (6) 

[03, a±] = ±a±. (7) 
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The abstract Lie-algebraic generalization of ([T]) is thus 14 1 



H 



+ {a{a- + a+) + l3a3)B + Hr,. 



(8) 



The simplest example of a representation of ([8]) is a single spin-1/2 particle interacting with a single one-dimensional 
harmonic oscillator (both su(2) and ho(l) are then given by irreducible representations). The Tavis-Cummings system 
is formally equivalent to several spin-1/2 systems interacting with a single harmonic oscillator (su(2) is then given 
by a reducible representation incorporating different spins, but ho(l) is still represented irreducibly) . The atom-field 
system we consider in the present paper is dual to the Tavis-Cummings model: its formal equivalent is a single 
spin-1/2 system interacting with several harmonic-oscillator wave packets. 

The reason for reducibility of the A^-atom representation of su(2) discussed in [l^ is mathematically very deep: 
multi-atomic Hilbert spaces are described by tensor products of single-atom representations, but tensor products of 
irreducible representations of su(2) are reducible. So, physically, it is the multi-particle structure of the A^-atom 
system that forces the representation of su(2) to be reducible. 

Now, for a quantum optician a cavity quantum field is an ensemble of many harmonic oscillators. The relevant 
algebraic structure is given by the multi-particle harmonic oscillator Lie-algebra ho(m), with m = oo not excluded, 



The standard representation employed in quantum optics has the form typical of irreducible representations, with 
ao(fc) = 1. But, in light of what we have written above, a multi-particle system should, in general, be represented 
reducibly. So, is ao(fc) = 1 really obvious? In what follows we will see that various regularizations that plague 
quantum field theoretic calculations may indicate something exactly opposite. In our opinion the issue is essential for 
a correct formulation of field quantization and is, in principle, testable in experiments with vacuum Rabi oscillations. 
Let us explain this viewpoint in more detail. 

There is a simple argument showing that if the harmonic oscillators forming the field are described by wave packets, 
the resulting representation of the harmonic-oscillator algebra is reducible. In order to show it, let us return to the 
Alicki-type Hamiltonian Hqj = a^a ® hoj (in fact, discussed already in p^). where uj = a;|a;)(tL'| is some operator, 
acts in the Hilbert space of states {tp) — ^il'n,u\'n',U!), spanned by the eigenvectors \n,Lu), 



Let us note that acts effectively as the usual Hamiltonian fnoa'^ a typical of the harmonic oscillator with frequency 
cij, but this is true only in the subspace spanned by \n,Lo), with fixed lo. An important difference between huja'^a and 
can be seen if one computes the average 



If \ip) is entangled, then different frequencies may be related to different numbers of excitations. Although describes 
a single harmonic oscillator, the average looks as if we considered an ensemble of oscillators with different frequencies. 
An interpretation of these facts is obvious: {ip) is a single-oscillator wave packet, and lo is not a parameter but an 
eigenvalue (a quantum number). So the Alicki-type Hamiltonian describes the usual quantum harmonic oscillator, 
but with quantized u. Quantization of uj becomes more natural if one recalls that typical us occurring in quantum 
harmonic oscillators are functions of observables (a magnetic field, or center-of-mass position operators, say). 
Now, let us rewrite Hc; as follows 



[a-_(fc),a+(fc')] 
[a±(fc),ao(fe')] 
[a3(fc),a±(fc')] 



<5fe,fe'ao(fc), 
[a3(fc),ao(fc')] =0, 



(9) 
(10) 
(11) 



Ht2;\n,iLj) = huja^ a\n,uj) = ntihj\n,uj). 



(12) 




(13) 




(14) 




(15) 



UJ 



where 



(16) 
(17) 
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The decomposition (fT5)) together with the fact that the eigenvalues of are nhw shows that Hq^ has many (but not 
all) properties typical of a quantum-field Hamiltonian. The operators a_(u;) — a^, a+{u!) = alj, a-sioj) = al^cii^, and 
ao{w) = 1 (8) form a representation of the harmonic-oscillator Lie algebra 

[a-{Lu),a+{uj')] = 6^^^'ao{Lj), (18) 
[a±{Lu),ao{io')] = [a3iio),aQ{u')]=0, (19) 
[a3(w),a±(w')] = ±S^,uj'a±{uj), (20) 

which is reducible since aQ^uj) is a projector and not a multiple of the identity. 

We have shown that a single-oscillator wave packet leads to a reducible representation of ho(m), where m is the 
dimension of the space spanned by the eigenvectors of uj. Let us make one point very clear already here: In this 
representation, m in ho(m) is not the number of oscillators, but the number of eigenvalues of uj in the one-oscillator 
wave packet. In the next Section we will show that a representation of ho(oo) corresponding to N single-oscillator 
wave packets, each of the oscillators existing in quantum superpositions of all the possible frequencies, is a natural 
candidate for an algebra of electromagnetic field operators. 

Before we launch on detailed calculations, let us briefly explain why 1 < iV < oo reducible representations of 
the harmonic-oscillator Lie algebra may be precisely what is needed for a well defined quantum field theory. Let us 
consider a single-oscillator Hamiltonian, but now with the vacuum-energy term included, i.e. 

= (21) 



J2f^{ala^ + lao{^))- (22) 



In such representations vacuum is represented by the entire Hilbert subspace of all the states that are annihilated by 
all Ooj. If the oscillators are bosons, the role of vacuum is played by a zero-temperature Bose-Einstein condensate, i.e. 
any state of the form 



|0^) = |0^)(8---® |0^), 


(23) 
(24) 


auj\Oip) = for all uj. 
:s has the Hamiltonian 


(25) 


N 


(26) 



where H^'^ is the Alicki-type Hamiltonian acting in the Hilbert space of a jth wave packet. The average energy of 
the Bose-Einstein condensate (i.e. the energy of a vacuum) is 

UJ UJ 

(27) 

where Z — maXi^H-^o.i^P}, and Xuj = \ipo.uj\'^/Z, lim^j^oo Xi^^ — 0, is the natural ultraviolet cut-off function resulting 
from normalizability, |'/'o,i^P = 1, of the condensate wave function. 

The cut-off appears only because the vacuum state is normalizable — it is hard to imagine in quantum theory 
a reason more fundamental. However, there is a price for it: vacuum states are in such theories non- unique (as all 
Bose-Einstein condensates) . Since the existence of a unique vacuum state is one of the Wightman axioms [3, 113] , 
quantum field theories based on this type of reducible representations of ho(oo) have to be non-Wightmanian. In this 
context it is worthy of mentioning that the non-Wightmanian aspects of the theory are not in conflict with Poincare 
covariance and gauge invariance of electrodynamics, but concrete technical forms of these conditions are quite different 
from what we are accustomed to (see [l9l|). 

In formula ([TT]) we encounter the characteristic product of two parameters, = ZN. This single parameter will be, 
effectively, the only free element in our reducible-representation treatment of Rabi oscillations. Let us note that the 
(thermodynamic) limit N oo with = const, implies shifting the cut-off to infinity, since Z ^ is equivalent to 
Zi^ —f for all ws, with Z^^ = 1. We will later see that although the thermodynamic limit makes vacuum energy 
divergent, the limiting form of vacuum Rabi oscillation is well defined and carries information about the value of 
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III. REDUCIBLE REPRESENTATION AND ITS DECOMPOSITION INTO IRREDUCIBLE 

COMPONENTS 



Let us assume that the electromagnetic field is a gas consisting of N indefinite-frequency bosonic oscillators of the 
type described in the Introduction. The basis in the iV-oscillator Hilbert space H is given by tensor products 

\nu,i . . .n^f,) = \ni,u!i) . . .\nN,U!N) (28) 

where the ws belong to the set of all the frequencies allowed by the cavity boundary conditions. The vacuum at zero 
temperature is assumed to be the pure state 

\0) = \0)...\0)= O^i---O.«|0c.i...0^„), (29) 

UJl . . .UJj\j 

\0) = Y^OM. (30) 
ElO-l' = 1- (31) 

UJ 

Taking a^^ in the form (jl6p , we assume that the atom-light system interacting with the reservoir will be described by 
the A^- wave-packet reducible representation of the Hamiltonian ([8|), where 



















«+ 






as 













®I®...®I+... + I®...®I®a^'j =a„, (32) 

®I®...®I+... + I®...®I®a}}j =aj,, (33) 

® . ..® I + . . . + I ® .. .® I (^ala^ = N^, (34) 

^(^{l <S) «)/«). ..«)/H h/(g)...«)/«)(l«) = I^, (35) 

is a reducible representation of the Lie algebra p ^ -([20 |) . and J±, J3 are given, as in by the spin-1/2 representation 
of su(2). The physical meaning of this representation can be inferred from the form of 03, which is the sum of number 
operators of N independent harmonic-oscillator wave packets, and for A'^ = 1 it reduces to the single-oscillator wave- 
packet representation from Section IL The representation was introduced in the context of quantum optics in p^ . 
Preliminary results on its implications for Rabi oscillations can be found in 
Let us denote by 'Huji...ujn subspace spanned by jn^^^ . . . n^j^). Obviously, 

H = (36) 

OJi , . ...ujn 

For any sequence of frequencies ujj^ and for any oj, one finds that the subspace 'Huji...ljn invariant under the action 
of ([32)1 - ([35)) ■ The central elements ao — satisfy 

I^\n^,...nuj^) = J^\n^^■■■n^^), (37) 

where s is the number of occurrences of lu in the sequence wj^. . . ,lun- As we can see, is the operator of frequency 

of successes, known from quantum laws of large numbers [2ll . [2^ . [23l . [23 |. Moreover, if s = 0, i.e. when none of 
(jji, . . . , uipf equals LU, then 

ajn^^ . . .n^j^) = o^ln^i . . . n<^„) = 0. (38) 

Note that there exist nontrivial states annihilated by creation operators. In consequence, for any -I- 1 different 
frequencies wi . . . Ci^Ar+i, one finds 

al ...n\ H^a,, a}, ...a}, H = • • • = a, , . . . a, , H = 0. (39) 

—UJl ^i)jv + l —oJi— 0)2 —UJn + 1 —0)1 —UJn + 1 \ ' 

The full system-reservoir Hilbert space Hs+r also can be split into a direct sum of subspaces invariant under the 
action of 

H = fi(|^T3+^i^ + g('7-at,+a+a^)) 

+ Ha^ + al)+PIl.)B + HR. (40) 
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Indeed, let Ha denote the 2-dimensional Hilbert space of atomic states, spanned by the excited and ground states, 
|e) and \g). Then 

Ws+fl = nA®n^,...u^®nR (41) 

where 

HnA®nu.,...^^®nR c nA®n^,...u^®nR. (42) 

The above observations, together with the known statistical properties of the operator of frequency of successes, 
supplemented by the fact that the right side of 

= (43) 

reduces in 'Ha®T^lji_...u:n ®^-R to multiplication by an appropriate s/N, make calculations in our reducible represen- 
tation as easy as those in the irreducible one. 

It it is interesting that the dynamics in each of the subspaces Ha ® "Huji... uj,^ ® H-r involves at most N different 
frequencies even if the number of frequencies admitted by boundary conditions is infinite. There are, of course, 
infinitely many such subspaces, corresponding to all the possible distributions of ujj into A'^ locations. 

In [l| we assumed the dynamics starting with the initial condition |e, 0)(e,0|. Its reducible analogue is 

p(0) = |e,0)(e,0| 

= |e)(e|® E 0^.---0.«0:, ...O:^|0^,...0^„)(0^;...0^^|. (44) 

The indexing sequence in 'Huji...uin uniquely determines the parameter s. Assuming temperature T = and (|44p as 
the initial condition, the resulting master equation (derived later) will involve only a subspace of Ti.uji...uK, namely, 
the one spanned by the following s + 2 orthonormal vectors 



|e,0^,...^„) = |e)|0„, ...0„...0^...0^„) = (45a) 

1^' ^LOi ...UJj^ / 



.lLV...c.„> = |g)|0.,...1....0....0.„) = |2), (45b) 



\9Ali...J = |g)|0„,...0„...1^...0^„) = |s + l), (45c) 
\g,Ou.^...u.^) - |g)|0^, ...0^...0^...0^„) = |s + 2). (45d) 

Clearly, \g, * € {1, • . . , s} denotes the state of the atom-field system in which the atom is in the lower state 

\g), and out of s oscillators whose frequency is tj, it is the ith oscillator which is excited to the first energy level. 



IV. DRESSED STATES IN H/ 



Denote by niji...t^jv the projector on Ha S?" Huji...ujk, and by tTuji...uik the one on the (s 4- 2)-dimensional subspace 
spanned by (|45|) . The total Hamiltonian consists of three terms: the Jaynes-Cummings atom-field part Hjc — fi^, 
the reservoir Hamiltonian Hr, and the system-reservoir interaction. The dressed states relevant for our problem are 
the eigenstates of 



(46) 



It is instructive to write in the basis (|45)) the matrix fl{s)ki ~ (fc|f2i^j...(j„ 



n{s) 



/ 2 Vlv \/lv ■ ■ ■ Vlv \ 
^ ^2 ... 
^0 I ... 



Vlv 
\ 







(47) 
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As we can see, ri(s) does not explicitly depend on the concrete values of frequencies that index ^uji...ujn^ but only on 
the number s of times the resonant frequency lo occurs in the indexing sequence uji . . . wjv The eigenvalues are 



UJ s 



no{s) = - 

with the corresponding orthonornial eigenvectors 



UJ 

'r 











Vi/ 



\n,{s)) = 



/ 1 \ 







V / 



\nk{s)) 



(48a) 
(48b) 
(48c) 



(49a) 



V J 



(49b) 



(49c) 



and 



|f7_(,s)) = 



/2s 



f-f\ 



V / 



, |f^+(s)) = 



/2s 



V J 



(49d) 



To avoid possible confusion let us remind that equations (|49p define an orthonornial basis in the relevant subspace 
of Ha ^ W[^i...wjvi and not in the full Hilbcrt space H. The dressed states are related to the bare ones ps)) by 



l^o(s)) = l5,0c^i...c^„), 



\n+is)) 
\n_is)) 

l^kis)) 



7f(7iEl5'iS...c.«) + le,o.,....„)), 



1 / 1 



7f(7iEi5-iS.....J-l^.o, 



(fc+l) \. L _ 1 
UJjv/' ft- — 1, 



(50a) 
(50b) 

(50c) 
(50d) 
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V. MASTER EQUATION AT ZERO TEMPERATURE 

We are looking for a density matrix p{t), acting in Ha 'S) "H, satisfying the initial condition p4)) . Since our goal is 
to compute the atomic ground-state probability, 

Pgit) = Tr\g)(g\pit) 

= J2 Tr|g)(.g|n^,...^„p(i)n„j...^^ 

UJ1...UJN Uj[...Uj'j^ 

= Tr|g)(5|n 

(51) 

we have to derive and solve an effective master equation for the projected density matrix 

Pu)1...u)n{^) ~ '^ljJl...l^JN P{^)^i^l---i^N ■ (52) 

The Jaynes-Cummings Hamiltonian can be analogously split into 

Hjc — E n^^i-.-WN-^ (53) 
= J2 (54) 

Each block H^-^^,,,^,^ can be yet further decomposed, 

Hui...uN = ^wi...wjv ^ (55) 

-^LJ1...WN ^ ""wi ...Wjv^t^l...lJN ■ (56) 

Let 7'(euji...aj„) denote the spectral projectors of 

-f^Wl-.-OJ^v — ^ ^ ^ ^OJi . . .CtJjY (^OJi . . .CtJjv ) ■ (57) 

The spectrum of H is inifinitcly degenerated since each H'^^ has eigenvalues described in the previous section, 
and there are inifinitely many sequences loi . . . luj^ corresponding to the same s. So let 

^(e) = E ^(e-i.-«)- (58) 

Taking into account the system-reservoir interaction Hamiltonian 

AB^{aia^+ql) + PNjB, (59) 
we find that the full density matrix p{t) satisfies at T = K the usual Markovian master equation [2^ 

p = -^[n,p]+Y,J{u;) (A{uj)pA\i.)-^[AHu;)A{cj),p]\ (60) 

where 

Aiu;) = E 'Pie)AVie'). (61) 

e' —e—u 

Since 

[7ra.i...^«:^] = [7r^i...i^«,'P(e)] = , ^] = (62) 
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we arrive at 



(63) 

with 

= '^-i...^„7'(e^,...^„)A7'(eL,...^„). (64) 

There are no jumps between energy eigenstates belonging to subspaces indexed by different sequences uji...u)m- 
Employing (jSOp and the explicit form of A, we obtain 

^sl^ms)){^+{s% (65a) 
^\n^{s)){n+{s)\ (65b) 

^^\n,{s)){^l^{s)l (65c) 

A^i...c.„(^^fc(s)-f^-(s)) 
= A^i...c.«(^^fc(s)~f^o(s)) =0, /c=:=l,...,s-l. (65d) 

An important and rather unexpected result is that jumps involving |r2fc(s)), fc = 1, . . . ,s — 1, are not allowed [26j . 
This is why the master equation for Puj^...ujtf{t) reads 







-f^o(s)) 


A-uJi.. 






A-uJi. 













2±_ 

N 



X Q|f7o(s))(l]+(,s)|p^,...^„|17+(,s))(17o(s)| - h\n+{s)){n+{s)\,p^,...^,] + 
+ 7(f7_(,s)-f^o(s))a2A 

X Q|0o(s))(f^-(s)|p.,....„|i^-(s))(i7o(s)| - \[\n^{s)){n,{s)\,p^,...^,u 

X Q|f7_(s))(r!+(s)|p.,....„|f7+(s))(l]_(s)| - h\n+{s)){il+{s)\,p^,...^,]+^ , (66) 

and, up to the presence of s/N in terms involving , is identical to the last equation from the Appendix of 
In consequence, we can directly apply the results from [l[ to (|66|) . One of the consequences of (|65)) is the time- 
independence of 

Plji...ujn P^l ■ ■ -t^Af (^) TrTT^^ __^jyp(^) 

= Trn^,...^„p(t) (67) 

Pui...u}N the probability of finding the sequence loi . . .uj^ if one randomly and independently selects each lu. The 
probability of finding uj equals = jO^^P, so Pi^i...cjjv = ^'^i ■ • ■ ^t^N- Now, let 

(68) 

be a normalized solution of (j66p . with the initial condition 

^ V+(s)>(^^+(5)l + |r!-(.s))(i^-(s)| - \n+{s))(n^is)\ - |i]_(s))(r!+(s)|^ 

(69) 
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The parameters 71 = j{n+{s) — rio(s))a^, 72 = j{fl-{s) — f2o(s))a^, and 73 = 7(ri+(s) — f2_(s))/3^/2 are related to 
the system-reservoir interaction Hamiltonian in a way that is identical to what was found in [l|. In order to have a 
well defined limit N ^ 00 we assume they are independent of s. However, in ([66]) 71 and 72 are additionally multiplied 
by s/N, a fact that introduces an s-dependence into 71 (s) = "fis/N, 72(5) — ^2s/N, keeping 73 independent of s. 
Applying the damping-basis method (see the Appendix), we find 

. , 1 1 71 (s) + 73 ^ 

Quii...un\}') — ~7i 7~\ I ^ ^ 

2 7i(s) - 72(5) + 73 



(71 (5) - 72 (s) + 73) 1^!+ (s)> (s) I + 73 (s)) (s) I 

+ (7i(5)-72(s))|r!o(s)>(r!o(,s)| 
l7i(.s)- 72(.) + 273„_^.^,^_(^)^^^_(^),_,^^^^)^^^^(^^, 



2 71 (s) - 72(3) -1-73 

+ |r!o(s))(r!o(s)| 

_ i^-2.,7I7]7t^- -^'-'+7'°'"-- *|f^^(,))(i74.)| 

- ie-^^/^*e--^^^^^'|^7_(.))(^^+(,s)|. (70) 

The conditional probability of finding the atom in its ground state, under the condition that the sequence is cji . . . cj^r, 
reads 



Pg{s,t) = Tr\g){g\g^,,,,^j^{t) 

^ ^ 1 7l(s) - 72(5) + 273 _22^t 1 7l(s)-72(s) _2iW+23t 

4 71 (s) - 72(5) + 73 4 71 (s) - 72(5) -I- 73 

--e J * cos2qy^s/Nt. (71) 

It is instructive to confront the formula (|71[) with the one we would have obtained in the formalism from [1] if, at 
T = 0, we employed an irreducible representation with ag — Zl, for some constant Z > 0, i.e. with [a, a^] = Zl (if 
Z < 0, then a is a creation operator). We first define a = aj\fz^ a' = /\fZ and then perform calculations with H 
ex pre ssed in terms of these new operators and appropriately rescaled parameters (see the Appendix) . The probability 
is [23] 

ij. 1 7iZ - 72Z -I- 273 _22£( 1 71Z-72Z _Ti£±23j 

" 4 7iZ - 72Z -I- 73 471Z-72Z + 73 

1 71^ + 72^ + 73 t „ , , 

--e 3 *cos2qVZt. (72) 

Obviously, Z occurs in (jTS]) in the same place as s/iV in ((7T|l . This is consistent with the fact that in subspaces 
characterized by s, the right-hand side of an analogous commutator involves s/N. All irreducible representations 
imply the same physical result provided one defines observable parameters by their renormalized forms: qph ~ Q\^, 

71, ph = 71-2^: 72, ph = 72-2, 73,ph = 73- 

Returning to the reducible representation, the ground-state probability we are looking for is the weighted sum 



PaW = P'^i-'^ivTr|g)(g|e^,...^„(t) 



N 



Ui(l-^.)^-^P.(s,i). (73) 

This is basically the final formula that should be compared with experiment. Before we do that, however, we have to 
relate the bare parameters 71, 72, 73, and q, to their physical, renormalized counterparts. 

In order to do so, we consider the asymptotic limit N 00. The weak law of large numbers (Feller's theorem 
dllli]) implies 

lim p (t) = 1 ^ ^1^" "^2^" + ^^3 ^-^* ^ liZ^~l2Z^ ^_ 7iz^+73 t 
N^oa ^ 4 7iZc^ - 72 2'^ -I- 73 471^0; - 72^i.j -I- 73 

1 71 +72 +73 * ^ I ^ 

-^e 2 Uos2q^/Z^t. (74) 
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The result is practically identical to (|72|) . q\fZZ is the effective physical coupling, where = \OJ\^ is the probability 
of finding the mode lo. Defining Z = maxi^jZi^}, we first of all note that the formula involves an automatic cut-off 
Xi^ = Z^^jZ^ < < 1- The observable parameters are identified with = gV^, 7i,ph = 7i^, 72, ph ~ I2Z, and 
73, ph = 73- The reducible-representation asymptotic formula reconstructs exactly the standard one following from the 
analysis given in fl] if we assume that Xo; = 1 for frequencies belonging to the optical regime. The formula for Pg{t), 
valid for all values of 1 < iV < 00, finally becomes 



s=0 



^ , _ 1 Wz (71, ph " 72, ph) + 273 ^_ T2,ph-ivrz ^ ^ 1 ^rz (7i,ph - 72,ph) ^_ ^i,ph:p-+T-3 j 
4 1^ (7i,ph - 72,ph) + 73 4 (7i,ph - 72,ph) + 73 



1 ife ("'l,ph+^2,ph)+-'3 



(75) 



The limit \\mM^ooPg{t) reconstructs (f72|) with Z = Z. The larger A'', the less important the exact value of Z. For 
N of the order of 10^ or higher, plots of Pg{t) are insensitive to changes of Z if the product ^ = NZ is kept constant. 
Such a limit, N — )■ 00, <r =const, is precisely a thermodynamic limit with fixed effective number <; of oscillators that 
interact with the two-level system [compare Eq. (|89p]. The thermodynamic limit implicitly removes the cut-off since 
Z — (^/N ~f 0, Z ~ maXi^{Z^}, Z^^ — 1, implies shifting cut-off to infinity. Requirements of mathematical 
consistency imply that N is finite. For physical reasons, however, N must be very large, and thus Z is small but 
non-zero. 

Let us now assume that the cavity is identical to the one employed in The mode has a Gaussian structure, so 
we have to correct ([75]) in a way described in detail in jjj . Denoting by d and w the cavity length and the Gaussian 
width, respectively, we get the effective probability 



N / 



s=0 



1 _ 1 ih (7i,ph - 72,ph) + 273 ^_ -'i-^^TT^ t _ 1 1^ (7i,ph - 72,ph) t 

4 Wz (7i,ph - 72,ph) +73 ijfz (7i,ph - 72,ph) + 73 

1 lfe(^l,ph+-'2,ph)+-'3 ^ / g ~| 

^^e- ^v/-™/^ cos2<Zp,^— 1|, 

(76) 

where t in (I76p is the effective time [l[. Analogously, the Gaussian- mode correction to the irreducible case reads 

~irr/,N 1 I71Z-72Z + 273 - -p.^.. t 1 71Z-72Z - :i^+7 t 

Pg [t) = 1 e 2v^u,/d g 2v^™/d 

4 7iZ - 72Z + 73 4 "fiZ - 72Z + 73 

__e 4v^™/d 'cos2qVZt. (77) 

Fig. 1 shows the time-dependence oi Pg{t) for iV = 10^ and various values of Even more suggestive is the plot of 
the difference Ip^it) ~ ^5(^^)1 (Fig- 2) compared with the error bars taken from the data of Brune et al. [2|. It is 
clear that this concrete experiment cannot discriminate between the limit iV — > 00 (i.e. the standard theory based on 
irreducible representations) and the alternative non-Wightmanian theory with any finite <; > 400. 



VI. ENERGY DECAY 



Identification of physical parameters with the renormalized ones is supported by the analysis of energy losses. In 
irreducible representations with [a, a'''] = Zl, the average energy of the atom- field system inside of the cavity, at 
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0.00002 0.00004 0.00006 0.00008 

t [s] 

FIG. 1: Probability Pg{t), Eq. ([76}, of finding the atom in the lower state \g) for ? = 1 (dotted), ? = 100 (dashed), and ? = 400 
(solid), with A'^ = lO'' for all the three cases, t is the effective time. The filled circles represent an analogous probability 
obtained for irreducible representations. Error bars taken from the experiment of Brune et al. provide a natural measure of 
distance between predictions of the alternative theories. For higher values of ?, say ? = 1000, the reducible representation 
becomes indistinguishable from the irreducible one. The remaining parameters are: q = 477rl0'' Hz, 71 = 72 = 83.912 Hz and 
73 = 0.07«7. 



T = 0, E{t) ~ hVl{ty'^'^ ~ hTTflp{t) [see the Appendix for the exphcit form of p(t)], is given by 



1 f 7iZ - 72Z /— 71Z - 72.Z + 273 1 _2l£±22.( 

qVZ — — > e 2 



2[ 71Z-72Z + 73 71Z-72Z + 73 

1 7iZ - 72Z + 273 / _Ta£4 

+ - — [uj-qVZje 2 . (78) 

2 7iZ - 72Z + 73 ^ ^ 

In [l[ we showed that best fits to experimental data are found for 71 =72. Inserting 7 = 71 = 72 into (|78p . we get 



n{tY" = -- + uje-^* + qVZe-^^e-^' ~ I). (79) 



It is evident that the energy damping parameter is jZ, and not just 7. 
The reducible-representation result is similar, 



If Tils) -72(5) , /^7l(s) - 72(5) + 273 1 _2ii£l±23t 

2 1 7i(s) - 72(5) + 73 V ^ 7i(s) - 72(5) + 73 J 
_^l7i(^) -72(.) +273 / _ fV\ ^.-j^, (30) 



2 71 (s) - 72(3) + 73 V V 
VL{s,t). (81) 
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FIG. 2: The difference |Pg' (t) —pg{t)\, with Pg^{t) and Pg{i) given by (|77p and (|76p . respectively, f is the effective time. ? = 400 
(sohd), ? — 1000 (dashed), ? = 5000 (dotted), A'^ = 10^. Stars represent the error bars taken from the experiment of Brune et 
al. The curves remain practically unchanged for higher A'', so the plots survive the thermodynamic limit with ? = const. 




FIG. 3: Comparison of r2(t)'" (dashed) and Q.{t) (solid), for (A) ? = 10 and (B) ? = 400 (curves are indistinguishable). The 
other parameters are: = 10^, 71, ph = 72, ph = O.lgph, 73 = O.OOlgph, gph = 477rl0^ Hz. 



As before, the right-hand side depends on s and not on the exact form of the sequence loi. ..lon- Repeatmg the 
reasoning from the previous sections, we find 

W) ^Y.{] ^'(1 - Zf^'WJ)- (82) 
Fig. 3 compares the two expressions for various values of ^, after having renormalized the parameters. 
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VII. VACUUM COLLAPSES AND REVIVALS FOR ? < oo 



The discussed Rabi-oscillation data do not distinguish between the two alternative forms of field quantization: one 
can always increase the value of ^ and produce a theory indistinguishable from the standard one within some given 
error bars. However, monitoring the oscillation long enough we can determine the value of provided < oo. 

For a finite N, the eigenvalues s/N are distributed around the most probable value s/N « Z (i.e. s « resulting 
in distribution of Rabi frequencies even in exact vacuum. In consequence, instead of a single-frequency oscillation we 
obtain beats analogous to those occurring in Rabi oscillations in the presence of a coherent light. Rabi frequency in 
exact resonance and in a subspace characterized by s equals 



2q^ - 2q^/Z^/^ = 2gph,/-, (83) 



NZ V <^ 



so that the most probable Rabi frequency is 2(7ph. The first revival of a collapsed vacuum Rabi oscillation will have its 
maximum when phases of neighboring and dominating oscillating terms differ by the factor of 2tt. This means that 
the revival time, tr, can be determined from 



2gphir - 2qpi,J- -U = 27r. (84) 



We find 



tr = (<?+ V^(<;-l))7r/gph. (85) 



If we take dissipation into account, the revival can be seen only if the oscillation occurring in ((75|) is still visible. The 
amplitude of oscillation is described, at s « ^, by 

1 Tl.ph+T2.ph + ^3 J 

£=-e-- ^ *^ (86) 



7 = 71, ph + 72, ph + 73 

= -- ln2e (87) 

t r 

with < £ < 1/2. Actually, due to interference effects the amplitude of the revival is some ten times smaller than e. 
The dependence of 7 on e is shown in Fig. [H Fig. [5] shows the revival of the decayed Pg{t). The parameters used in 
Fig. [1] imply e ~ 10~^°, so the effect would not be visible in experiments where error bars grow with time similarly 
to those from Moreover, the data monitor oscillations for four Rabi periods, whereas for — 400, assuming the 
most optimistic scenario, the first revival should be seen after approximately 800TRabi. 



VIII. FINITE N OH N ^00? 



Cavity boundary conditions imply that there are infinitely many ujs. The probability of finding a given to is 
Zoj = |0(jp- Physical intuition suggests that Z — maxi^{|Oi^p} should be a small but non-zero number. Summability 
of probabilities to 1 implies that Z^ — > if w — *■ 00. On the other hand, evidently Zi^ « Z for optical frequencies since 
Xu! = Z^^/Z plays a role of a cut-off which, in the optical regime, should satisfy Xt^ ~ 1 (see, for example, a discussion 
of this point in p^). 

One of the important properties of field quantization in terms of reducible representations is the correspondence 
principle with standard regularized quantum optics (it is beyond the scope of the present paper, but computation of, 
say, resonance fluorescence in the reducible-representation formalism indeed introduces in those places where one 
puts cut-off by hand in the standard formalism). What is interesting, the role of correspondence principle is played 
by the weak law of large numbers. 

The weak law, iV — 00 with Z^ — const, replaces all s/N by probabilities Z^. At the level of representation, the 
weak law follows from spectral representations of central elements aQ{Lo) — occurring at right-hand sides of the 
commutator [g^,g^,] = S^^'L^j (recall that are the frequency-of-success operators known from quantum laws of large 
numbers). The reducible representation may be regarded as a "quantized" form of the standard naively regularized 
irreducible representation [a^, al^,] = Sujui'ZljI, where Z^ is a regularizing function. The naive regularization is known 
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FIG. 5: Probability of the atomic ground state as a function of the effective time t for 71 = 72 = 83.912 Hz, 73 — 10 Hz, and 
? = 400 (these parameters yield e ~ 0.23). The revival time (|85|l . tr ~ 0.017 s, is indicated by the vertical line. It approximately 
determines the moment of maximal visibility of the revival. 



to be in conflict with Poincare covariance of canonical commutation relations. The replacement of the function 
by the operator leads to the correct behavior of the commutator under Poincare transformations jH, [IS, [sH, . 
The thermodynamic limit, TV — > 00 with ^ — NZ — const is different. Since Z = <j/iV, the bare charge satisfies 

qph [N 



The Jaynes-Cummings interaction [see ([31])] 

qa+a^ + H.c. = \j^qph(^+ {a^ (g) I (g) . . . (g) I + . . . + I (g) . . . (g) I (g) a^] + H.c. 

shows that the physical coupling parameter is, in fact, independent of A^: 

gph ^ _Q_ 



(^a^g) I g) . . .<Si I + . . . + I g) . . .g) I ^aj^ + H.c. (89) 



(90) 
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for any N. Let us stress, that although we do not yet have exact results on the existence of the thermodynamic limit, 
all numerical experiments show quick convergence, with growing N, of Pg(t) to a function whose shape is characterized 
by i;. 

Thermodynamic limits typically make some quantities divergent whereas some other quantities are well defined 
(think of a glass of water — the density of water will correctly behave in the limit, but the water mass will become 
infinite). A similar situation is encountered in our treatment of cavity QED. The probability Pg{t) becomes well 
defined even if N tends to infinity with kept constant, but the vacuum energy diverges. So, physically we have to 
assume that N is large but finite, similarly to Z which is small but nonzero. 

IX. OVERVIEW AND CONCLUSIONS 

All quantum optical experiments suggest that the idea of replacing classically oscillating amplitudes, ai^e"*'^*, by 
quantum harmonic-oscillator operators, Oi^e"*'^*, is correct. What is important, the representations of harmonic- 
oscillator Lie algebras we know from quantum mechanics textbooks correspond physically to quantum oscillators 
whose w is a classical parameter. (Heisenberg, Born, and Jordan were not aware, in 1925, of the role of eigenvalues 
for quantum theories — the point was understood by Schrodinger in 1926 [s^].) Simultaneously, for physical reasons 
the frequency uj is in actual oscillators a function of observables (for example, endpoints of a nano-scale pendulum are 
given by atomic center-of-mass wavepackets and, hence, its length is unsharp — the pendulum exists in a quantum 
superposition of different lengths I or, equivalently, of the corresponding frequencies lu = \J gjl). 

Attempts of taking this observation into account lead to oscillator wave packets existing in superpositions of different 
ws. The frequencies are given by eigenvalues of an operator a) that, in simplest models, commutes with other oscillator 
observables (in yet more realistic cases the operators tD will not commute with other observables). Since all ujs can be 
associated with a single oscillator, there is no relation between the number (typically infinite) of field modes and the 
number N of oscillators (that can be finite). 

For those who are trained on the usual Heisenberg 1925 oscillator, the conclusion seems counterintuitive: Now a 
single oscillator may contain many modes, and a single-mode problem may involve many oscillators. An algebraic 
implication is that once we replace w by cD, a new representation of the harmonic oscillator Lie algebra arises. The 
representation is reducible, and the commutator [a(j,ajj] is no longer proportional to the identity, but is a nontrivial 
operator that commutes with all the operators of the algebra. The vacuum average of this operator plays a role of a 
regularizing function that, at least in all the applications considered so far, appears in exactly those places where one 
expects regularization to occur. Simultaneously, at the level of spectra, instead of regularizing functions one finds 
the eigenvalues s/TV, s = 0, 1 ... TV, of the operator [a^^, a\\. The eigenvalues s/7V are not distributed uniformly. The 
weak law of large numbers implies that for a large N they concentrate around the probability Z^j — \OJ^ of finding 
the eigenvalue w in vacuum. The maximal value Z = max^lZi^} is a parameter that characterizes the structure of the 
vacuum state. The product <; = ZN represents an effective number of field oscillators that interact with the two-level 
atom. Since for physical reasons Z should be small, then keeping <; fixed we have to assume that N is large. Actually, 
so large that the thermodynamic limit, with <; = const, should be applicable. In relativistic version of the theory both 
Z and N are relativistically invariant. 

As we can see, there are reasons to regard the reducible representation of harmonic oscillator algebras as physically 
justified. Now we face two questions. First, can the available data really prove that the representation of the Lie 
algebra is the one we know from standard quantum optics? The answer is: Probably no, since the limit — > cx) is 
a correspondence principle with the standard theory, so one can always increase N to approach the standard theory 
within given experimental error bars. So the second question is: Does there exist a phenomenon whose observation 
could prove that N and <; are finite? The answer is, in principle, positive. For example, monitoring a decayed Rabi 
oscillation in exact vacuum we expect its revival after, approximately, 2^rRabi- Current data suggest ^ > 400. 

It is also worthy of mentioning that the parameters Z and N alone may not be measurable, since Z gets hidden 
in the renormalized charge Cph = ev^ (note the analogy between Z and the renormalization constant Z3). Once 
incorporated into renormalized parameters, Z disappears from large-A'^ asymptotic formulas, but what remains are 
the "cutoff functions = Z^/Z. The thermodynamic limit turns out to remove the cutoff. In effect, the standard 
steps of the renormalization procedure are here interpretable in physical terms. If we could confirm that physical 
and N are finite, we would have a proof that physical vacua are more similar to realistic Bose-Einstein condensates 
than to "unique and Poincare invariant cyclic vectors" typical of Wightmanian axiomatics. 
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Appendix A: T = solution in irreducible representations 



For T ~ 0.8 K the proportionality factor e linking transitions up and down between dressed states from the same 
manifold equals exp(— 2fi,q/fcT) « 0.999 997, but vanishes for T = (compare Fig. 10 in [l[). The case T = was 
not interesting from the point of view of the analysis given in [l|, so we made there the approximation e = 1. In 
consequence, an appropriate T — formula for p^g^{t) — Tr|(7)(g|yo(t) is missing in [l[. Below we give a detailed 
derivation of the solution corresponding to T = 0. This supplements [T], simultaneously showing how to solve the 
master equations occurring in the present paper, since mathematically they are identical. 

The irreducible representation is [a, a^] = 1; if [a, a^] = Zl, Z > 0, we redefine a/VZ a, a) j^fZ . The 
T = master equation reads 

p = -z[r!,p]+7i|i|f^o)(r!+|p|r!+)(r!o|-i[|r!+)(r!+|,p 

+72 |i|f^o)(f^-|p|f^-)(f^ol - 



+73 \ - 



(91) 



Eg. ([91]) . written as p = Lp^ defines the operator L. The eigenvectors Lpj — h^jPj define the so-called damping bases 

71 + 73 



Pi 

P2 
P3 
P4 



-(71 -12+ -f3m+){Q+\+j3\n_){n_\ + (71 - 72)|r!o)(f^ol, Ai 



(7172 + 7273)|^^0>(f^o|, A3 



0, 



= |f7+)(r2_|, A4 = 



71 + 72 + 73 



P5 = |f7+)(r2ol, A5 = -i(f7+-f7o)- 

P7 



r!_)(r!o|. As = -z(fi_-r!o)-^ 



71 + 73 



72 



A7 = - r2_) - 



71 + 72 + 73 



|r!o>(r!+|. As = 



P9 



71 + 73 



72 



The initial condition 



p{0) = \e,o){e,o\^^\n+){n+\ + ^\n^){n^\-^\n+){n_\-^\n_){n+\ 



1 



1 



1 7i - 72 + 273 
-Pi + o ; P2 



1 



1 



1 



2 71 - 72 + 73 2 71 - 72 + 73 



—7 : rP3 - xP4 - 7;P7 

72(71+73) 2 2 



(92a) 

(92b) 
(92c) 
(92d) 

(92e) 
(92f) 

(92g) 
(92h) 
(92i) 

(93) 
(94) 



leads, in exact resonance, to 
Pit) 



1 



1 



2 71 - 72 + 73 



Pi 



1 71 - 72 + 273 _22t 

-e 2 p2 



1 



2 71-72 + 73 



1 _2iot _2i±22±23* 1 2iat _2i±22±23f 

irr/,N _ 1 71 - 72 + 273 -42.t 1 71-72 _TL+23 

Pq (1^) — 1 . e . e 



72(71 + 73) 

j 1 f 



P3 



4 71-72+73 



4 71 - 72 + 73 



cos2gi. 



(95) 
(96) 



Appendix B: Quantum wealt law of large numbers 



Consider some projector = |fc)(fc| and a state vector = J2k^k\k)- The probability of observing a property 
represented by Pk is \ipk\'^ and the eigenvalues of are ("failure") and 1 ("success"). Now let us consider a system 
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of N copies of \ip), i.e. 

I*) - IV') <^ • • • «) IV') • (97) 
V ' 

N 

The operator that counts the number of successes in N independent measurements of Pk is 

Nk P/c«)l«)---«)H + l«)---«)l«)Pfe (98) 

Denoting Pi = P^, Pq = I — Pk, obtain spectral representation of in the form 

N 

Nk ^ H sPa,<»---<»Pa^, Aj^0,1- (99) 

s=0 AiH hAiv=s 

Let / : [0, 1] ^ i? be some function. Spectral theorem defines f{Nk/N) as 

N 

f{Nk/N) = J2 E /(^/^) ® • • • ® (100) 

s=0 AiH |-Aiv=s 

For N ^ oo the average 

(vI/|/(7Vfe/iV)|vI/) ^ E/(s/iV)(^)(|^fep)^(l-|V;,p)^~^ (101) 

converges to /(iV'feP) (by virtue of the law of large numbers for the binomial distribution). The latter implies also 
that 

hm II /(7V,/iV)|vI/)-/(|7^,f)|vI/) 11=0. (102) 

An analogous proof can be formulated for projectors Pk projecting on higher dimensional subspaces, for example for 
Pk = 1 |fc) (fc| . Attempts of formulating a quantum analogue of the strong law of large numbers lead to infinite tensor 
products, A'^ = oo, and thus to nonseparable Hilbert spaces. Field quantization in terms of reducible representations 
of the harmonic oscillator Lie algebra involves only a finite A^, even if the number of frequencies allowed by given 
boundary conditions is infinite. 



Appendix C: Naive regularization versus relativistic covariance 



Let a{k) be a scalar field annihilation operator that transforms covariantly under Lorentz transformations: 
U{Lya{k)U{L) = a{L-^k), with U{LyU{L) = 1. The standard commutator [a(fc), a(/c')1 = ZS{k - k')l, where 
Z is a constant and 6(k — k') a relativistically invariant Dirac delta, is relativistically covariant. The proof is trivial: 
On one hand, U{L)'^[a{k),a{k')^U[L) = [U{L)^ a{k)U{L),U{L)^ aik'^'^UiL)] = [a{L-^k),a{L-^k')^. On the other 
hand, U {L)'^ Z5{k — k')l U{L) — ZS{k — k')l — Z5(L~^k — L~^k')l since delta is invariant by assumption. The proof 
fails if Z is replaced by a nontrivial function Z{k), i.e. in naive regularizations. 

However, since U{L)'^\k) = \L^^k), a replacement of Z{k) by I{k) = 1 (g) |A:)(A;|, or by a similar object, removes 
the difficulty: U [L)"^ I {k)U {L) = I{L-^k). Simultaneously, in all averages the weak law of large numbers effectively 
replaces I{k) by Z{k) (see Appendix B). Another conclusion is that, as opposed to naive regularization, spectra of 
Hamiltonians do not depend on Z{k), but on relativistically invariant eigenvalues s/N (think of the dressed states we 
have described in the present paper). This is why, following Finkelstein [33|, we speak of regularization by quantization. 
For more details on relativistic aspects see e.g. |19| . 
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